Characteristic BEC scaling close to Quantum Critical Point in BaCuSi206 
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We report an experimental determination of the phase boundary approaching the quantum critical 
point separating a quantum paramagnetic state and the proposed spin Bose-Einstein condensate of 
triplons in the spin gap compound BaCuSi2 0s. The ordering temperature is related to the proximity 
to a quantum critical point at the lower critical magnetic field H c i = 23.52 ±0.03 T by a power law 
parameterized by critical exponent v. We obtain an experimental estimate of v = 0.63 ± 0.03 down 
to a temperature of 0.61K, which is in good agreement with the mean field prediction of v = | for 
the Bose-Einstein condensation universality class. 



Several spin gap compounds, including 
BaCuSi 2 6 [l,2], TlCuCl 3 [3-7] and Sr 2 Cu(B0 3 )2[8], 
have a singlet ground state in zero magnetic field with 
a gap to the lowest excited triplet state. The spin gap 
can be closed by an applied magnetic field, such that 
a quantum critical point (QCP) at a magnetic field 
H c ^ separates the quantum paramagnetic state from a 
state characterized by long range magnetic order. The 
order parameter for this transition is the creation 
operator for a triplet state. In the absence of U(l) 
symmetry-breaking anisotropy, the transition from the 
quantum paramagnetic state ((b*) = 0) to an ordered 
state with broken U(l) symmetry ( (6+ ) oc + iM^ 
where Mi and M^L are the x and y components of the 
staggered magnetization in the plane perpendicular to 
the applied field) may be interpreted as a Bose Einstein 
condensation (BEC) of triplons [9]. The magnetic field 
functions chemical potential, thereby providing a 

convenient means to tune a BEC to criticality at the 
QCP. 

The spin gap system BaCuSi2C"6 is unique because it 
provides experimental access to a quantum phase transi- 
tion of this nature in a spin lattice that can be described 
by a U(l) rotationally symmetric spin Hamiltonian [1]. A 
quantum phase transition that belongs to the same uni- 
versality class is observed in superfluid 4 He[10-12], but 
with a phase diagram that is comparatively inaccessible 
for experimentation [13]. 

The proximity to the QCP is expected to be related 
to the ordering temperature (Xc) by a power law Tq ~ 
(H — H c ]) u [14], which can be expressed in reduced form 



t = f(h) x (1 



where t 



T max 



hf (1) 
(-ffmax and Tmax 



represent the point on the phase boundary halfway be- 
tween H c i and i? c 2, the field at which the magnetization 
saturates) and f (h)h=i is finite. The mean field critical 
exponent v 



| is characteristic of the Bose-Einstein 



condensation universality class, and describes the scal- 
ing behavior of a 3D dilute interacting Bose gas near 
the QCP. The mean field estimate is appropriate since 
the upper critical space dimension of 2 (d c = 2 because 
z = 2 for this universality class) is exceeded [9, 15-17]. 

In this paper, we present a set of experiments on 
BaCuSi206 that examine the critical scaling of this phase 
transition in the vicinity of the quantum critical point. 
The experimental results are consistent with the v = | 
BEC critical exponent, and agree with Monte Carlo sim- 
ulations in the lowest experimentally accessible temper- 
ature window. 

BaCuSi2C>6 has a well characterized quasi 2D struc- 
ture [1,2,18] consisting of vertical Cu 2+ dimers arranged 
on a square bilayered lattice staggered between bilay- 
ers. The inter-layer, interdimer and interbilayer exchange 
couplings have been estimated from high field magne- 
tization data [1] to be J = 4.45meV, J' ~ 0.58meV, 
J" ~ 0.116meV respectively [1,2]. Representing each 
dimer by a pseudospin with two possible states (lowest 
energy states s z = singlet and s z = 1 triplet), and using 
appropriate transformations [1], the effective Hamilto- 
nian describes a 3D gas of hardcore bosons with in-plane 
nearest neighbor hopping and repulsive interactions. At 
low temperatures, the system undergoes a second order 
phase transition, which has been interpreted as triplon 
condensation [1]. 

Here, magnetic torque, magnetocaloric effect and spe- 
cific heat measurements are performed to obtain points 
on the phase boundary into the ordered state. Features 
in these thermodynamic quantities characterise the clas- 
sical 3D-XY phase transition into the ordered state at 
finite temperatures. 

Single crystal samples of BaCuSi2C>6 grown by a flux- 
growth technique [19] are used for these experiments, 
whereas previous measurements on this material [1,2,18] 
used single crystals grown by a floating zone technique. 
Flux grown crystals are chosen because of a lower impu- 
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rity content, a clearly defined Schottky anomaly in the 
zero field heat capacity, and narrower nuclear magnetic 
resonance lines. The characteristic features of the mag- 
netic ordering transition observed previously are identical 
in these samples, indicating sample independence. 

The specific heat of the sample is measured at 36T 
in a 4 He cryostat in the hybrid magnet at Tallahassee. 
The inset to figure la shows the characteristic lambda 
anomaly observed in specific heat, indicating a second or- 
der phase transition into the ordered state. The shape is 
identical to that in [1], which has been fit using directed- 
loop Monte Carlo simulations. The ordering temperature 
at 36T is plotted on the phase diagram in figure la. 

Magnetic torque measurements enable a sensitive 
probe of the magnetic ordering transition at low tem- 
peratures. They are performed in static magnetic fields 
up to 33T in a 3 He refrigerator in Tallahassee. Samples 
are mounted on the moving plate of a phosphor bronze 
capacitance cantilever, attached to a rigid plate rotat- 
able about an axis parallel to the axis of torque and per- 
pendicular to the applied magnetic field. The sample is 
mounted with a small angle (< 10°) between the applied 
field and the normal to the sample plane (easy axis c), 
such that the applied field exerts a torque on the crystal 
due to the difference in g-factor between the & and c ori- 
entations with g a = 2.053 ±0.007 and g c = 2.303 ±0.003 
[19]. The anisotropy in g results in an anisotropy in 
ff c j = (where A is the spin gap) [19]. Hence, on 
entering the magnetically ordered phase with increasing 
magnetic field, the anisotropy in magnetization causes a 
sudden increase in torque, as the field attempts to align 
the c axis more closely with the applied magnetic field. 
Torque measurements are made during magnetic field 
sweeps across the ordering transition at different tem- 
peratures. 

The signature of the ordering transition is seen in field 
dependent torque curves in a temperature range 0.61K 
- 3.3K (sample curves shown in figure lb). The field at 
which the phase transition occurs is obtained from the 
position of a sharp feature in the second derivative of the 
torque (example shown in the inset to figure lb) [20]. The 
feature becomes weaker at higher temperatures, but can 
be extracted up to T = 3.3K. Examples of the ordering 
transitions thus obtained are indicated by solid symbols 
on the torque curves in figure lb. Points on the phase 
diagram obtained from torque measurements are shown 
as solid circles in the phase diagram plotted in figure la. 

The magnetocaloric effect describes the temperature 
change of a magnetic material associated with an external 
magnetic field change in an adiabatic process. An abrupt 
change in the temperature with changing magnetic field 
indicates a large field variation of the isothermal mag- 
netic entropy, and is associated with an ordering transi- 
tion. These measurements are a good probe of ordering 
transitions in a rapidly changing magnetic field, so can be 
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Figure 1: (a) Points on the phase boundary are determined 
from magnetic torque (solid circles), magnetocaloric effect 
(open circles), specific heat (solid square) measurements and 
Monte Carlo simulation (+ symbols, dotted line is a guide to 
the eye). The inset shows the lambda anomaly in specific heat 
measured at a magnetic field of 36T. (b) Torque measured as 
a function of rising magnetic field at sample temperatures as 
indicated. The ordering transitions determined from a sharp 
feature in the second derivative (shown in the inset) are indi- 
cated on each of the torque curves, (c) Sample curve indicat- 
ing the temperature change due to the magnetocaloric effect 
measured as a function of up and down sweeps in a magnetic 
field. The ordering transition is indicated on the curve. 



used to trace the phase boundary [21]. Magnetocaloric 
effect measurements are carried out in fields up to 45T 
in a 4 He cryostat in the hybrid magnet in Tallahassee. 
Temperature changes are detected during magnetic field 
sweeps across the ordering transition at different temper- 
atures (shown in fig lb). The peak in lattice temperature 
at H c i {H c 2) in a rising (falling) field indicates a drop in 
magnetic entropy with ordering. Similarly, a dip in tem- 
perature in a falling (rising) field indicates the transition 
out of the ordered phase at H c ^ (H c 2). The position of 
the ordering transition for H < H max is obtained from 
the onset of the peak (as defined by the maximum in 
the first derivative) in a rising field and for H > H max 
from the onset of the dip in a falling field. The ordering 
transition thus obtained for a representative field sweep is 
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shown in figure lb. Points on the phase diagram obtained 
from magnetocaloric effect measurements are shown by 
open symbols in the phase diagram plotted in figure la. 

The Monte Carlo simulations for this system are per- 
formed using the directed-loop algorithm [22] (results are 
represented by + symbols in figure la). Estimates of 
interdimer exchange coupling are refined from [1], with 
revised values of J' = 0.51meV and J" = 0.168meV yield- 
ing better agreement of the Monte Carlo simulations with 
experimental points on the phase diagram. 
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Figure 2: (a) The circles represent estimates of obtained 
from fitting the lowest few experimental points on the phase 
boundary in figure lc in a window of increasing size £w, to 
eqn. 2 for different fixed values of v. The x axis t w labels the 
highest reduced temperature of the fit window. The dotted 
lines show the linear convergence of ' H' c \ values at t w = 0. 
The triangles represent estimates of ' H' similarly obtained 
from Monte Carlo simulation data for corresponding fixed val- 
ues of v, and similar convergence is observed (b) The circles 
represent estimates of v from fitting the lowest few experimen- 
tal points on the phase boundary in figure la in a window of 
increasing size tw, to eqn. 2 with H c ^ = 23.52T determined 
from figure 2a. The error bars are due to experimental uncer- 
tainty in determining values of He, leading to an uncertainty 
in the value of -ff c i . The triangles represent estimates of v 
from a similar fit to Monte Carlo simulation data. 

Experimental access to the particle-hole symmetric re- 
gion in BaCuSi206 enables us to extend the region near 
the QCP in which the power law can be fit. Eqn. 1 
describes scaling near H c ^ while t oc (1 + h) v describes 
scaling near i? c 2- In other words, the particle-hole sym- 



metry of the system implies that t is a function of h 2 : 

t = g(h 2 ) x [(1 - h)(l + h)] u = g(h 2 ) x (1 - h 2 Y(2) 

where g(h 2 ) varies more slowly than f(h) in the vicinity 
of the QCP. 

In this set of experiments, points on the particle-hole 
symmetric phase diagram first reported in [1] are refined 
using measurements on a flux-grown single crystal and 
extended to low temperatures close to the QCP (figure 
la). Critical scaling near the QCP is described by eqn. 
2. Importantly, we have access to sufficient data points 
in the low temperature region near H c ^ to obtain an 
estimate of the critical scaling exponent by fitting a power 
law dependence in that region. 

The power law dependence in the quantum critical re- 
gion is extremely sensitive to both the fit temperature 
range, and to the estimate of H c ^ [17]. We use an em- 
pirical convergence approach to determine the best esti- 
mate of -ff c i- Figure 2a shows the trend in the estimate 
of H c i (denoted as 'H' ^) obtained by fitting the low- 
est few experimental points on the phase boundary in 
figure la in a window of increasing size t w to eqn. 2 
(where g(h 2 ) is assumed to be constant to a first approx- 
imation) for different fixed values of v. Near the QCP, 
it is empirically observed from linear extrapolation to 
£ w = 0, that estimates of 'H' become less dependent on 
v, and converge to a single value irrespective of the value 
of v (figure 2a). Similar convergence to H c i is observed 
for the Monte Carlo simulation results (figure 2a). The 
convergence is due to the fact that the QCP is at -ff c x, 
independent of the path along which it is approached 
(characterized by v). From 'H' convergence, we obtain 
an estimate of H c i ~ 23.52T. The uncertainty in this 
value of H c i is estimated to be 0.03T, arising from the 
experimental uncertainty in measuring H c . The value 
of — 23.52 ± 0.03T thus obtained is then used to 
estimate the critical exponent v. 

The critical exponent v is estimated from fitting eqn. 2 
(with g{h 2 ) w const.) to the narrowest temperature 
range near the QCP with a statistically significant num- 
ber of experimental datapoints. Figure 2b shows the vari- 
ation in v with the size of the temperature window that 
is fit to eqn. 2 (points on the phase boundary are fit 
from the lowest value of t = 0.16 to a highest value of 
t = t w ). As the critical region near the QCP is ap- 
proached with narrowing window size, v approaches the 
theoretical mean field value of |. In the lowest experi- 
mentally accessible temperature window down to 0.61K, 
we obtain a value of v = 0.63 ± 0.03. Performing a sim- 
ilar analysis for data points from the Monte Carlo sim- 
ulation reveals the expected increase in v to the theo- 
retical mean field value as the temperature window is 
further reduced below currently accessible experimental 
temperatures (figure 2b). The experimental estimate of 
v = 0.63 ± 0.03 based on measurements down to tem- 
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perature t = 0.16 is consistent with the theoretical mean 
field prediction of v = | to within experimental error. 
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Figure 3: Points on the phase boundary determined from 
magnetic torque (solid circles), magnetocaloric effect (open 
circles) and heat capacity (solid square) measurements for 
H cl = 23.52 T, #max = 36.12 T, T m ax = 3.70 K and from 
Monte Carlo simulation (crosses). The lines represent eqn. 2 
with v = |. 

Figure 3 shows the comparison between the Monte 
Carlo simulation and the experimental data for H c i = 
23.52T. The lines represent the power law (eqn. 2) with 
v = |. The Monte Carlo simulation and experimental 
data are in good agreement in the lowest experimentally 
accessible temperature window. We empirically observe 
that the deviation from the power law at higher temper- 
atures is less in the case of the experimental data than 
for the Monte Carlo simulation. 

In summary, we performed magnetic torque and mag- 
netocaloric effect experiments to map out the phase dia- 
gram in the vicinity of the QCP in the spin gap system 
BaCuSi2C"6. Points down to 0.61K are fit to the power 
law t = g(h 2 ) x (1 - h 2 ) 1 ' (with g(h 2 ) w const.) to give 
a value v — 0.63 ± 0.03. The estimate we obtain for v 
near the QCP is close to the theoretical value v = |, 
whereas previous experimental measurements of v in the 
spin gap system TICUCI3 have resulted in lower values 
in the range 0.43 — 0.60 [3-7]. The other known exper- 
imental measurements of the BEC critical exponent v 
have been on 4 He adsorbed in aerogel [10-12], which is 
a realisation of a dilute Bose gas, but is experimentally 
limited by the presence of a random external potential. 
BaCuSi20e is a unique U(l) symmetric spin gap mate- 
rial that enables experimental access to a QCP separating 
a quantum paramagnet from a Bose-Einstein condensate 
[1] . It provides a novel experimental realization of a BEC 
in a grand canonical ensemble in the absence of an exter- 



nal potential, with the region around the QCP accessible 
by a tuneable external magnetic field. 
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